PROCEEDINGS 

™ SCIENCE 



Simulation of supersymmetric models on the lattice 
without a sign problem 



o 

< 



> 
O 

o 



> 

• I— I 

X 



David Baumgartner* and Urs Wenger* 

Albert Einstein Center for Fundamental Physics 

Institute for Theoretical Physics 

University of Bern 

Sidlerstrasse 5 

CH-3012 Bern 

Switzerland 

E-mail: baumgart@itp . unibe . ch, wenger@itp . unibe . ch 



Simulations of supersymmetric models on the lattice with (spontaneously) broken supersymmetry 
suffer from a fermion sign problem related to the vanishing of the Witten index. We propose a 
novel approach which solves this problem in low dimensions by formulating the path integral on 
the lattice in terms of fermion loops. The formulation is based on the exact hopping expansion of 
the fermionic action and allows the explicit decomposition of the partition function into bosonic 
and fermionic contributions. We devise a simulation algorithm which separately samples the 
fermionic and bosonic sectors, as well as the relative probabilities between them. The latter then 
allows a direct calculation of the Witten index and the corresponding Goldstino mode. Finally, we 
present results from simulations on the lattice for the spectrum and the Witten index for .A' = 2 
supersymmetric quantum mechanics. 



The XXVIII International Symposium on Lattice Field Theory 
June 14-19,2010 
Villasimius, Sardinia Italy 



* Speaker. 



© Copyright owned by the author{s) under the terms of the Creative Coinmons Attribution-NonCoinmercial-ShureAHke Licence. 



http: //pos. sissa. it/ 



Simulating supersymmetric models without a sign problem 



David Baumgartner and Urs Wenger 



1. Introduction 

Supersymmetry (SUSY) is thougiit to be a crucial ingredient in the unification of the Standard 
Model interactions as well as in the solution of the hierarchy problem. On the other hand, we 
know that the low energy physics is not supersymmetric, and consequently the SUSY must be 
broken at low energies, either explicitly or spontaneously. Since the origin and mechanism of 
spontaneous SUSY breaking is a non-perturbative effect, it can not be understood in perturbation 
theory - instead, non-perturbative methods are required. One way to study non-perturbative effects 
in quantum field theories is provided by the lattice regularisation. However, the lattice discretisation 
comes with various problems, for example the explicit breaking of Poincare symmetry, or the 
absence of Leibniz' rule, and it is therefore not clear at all whether and how SUSY can be reahsed 
within the lattice regularisation. 

If the lattice discretisation, for example, enjoys some exact symmetries that allow only irrel- 
evant symmetry breaking operators, which become unimportant in the infrared regime, so-called 
'accidental' symmetries may emerge from a non-symmetric lattice action and the full symmetry 
develops in the continuum limit. This is for example the case for (Euclidean) Poincare symmetry 
in lattice QCD or SUSY in = I SU{N) Super- Yang-Mills theory. In the latter case, the only 
relevant symmetry breaking operator is the gaugino mass term which violates the Z2N chiral sym- 
metry. Therefore, a chirally invariant lattice action forbids such a term and SUSY is automatically 
recovered in the continuum limit. For SUSY theories involving scalar fields, however, such a way 
out is not available: the scalar mass term m^|0|^ breaks SUSY, and there is no other symmetry 
available to forbid that term. In principle, even in such cases, some symmetries can be obtained 
in the continuum by fine tuning the theory with counterterms. The restoration of chiral symmetry 
for Wilson fermions is one such example. For SUSY, such an approach is in general not practical, 
but in lower dimensions, when theories are superrenormalisable, it sometimes is [1]. Yet another 
approach to SUSY on the lattice is to look for an exact lattice realisation of a subalgebra of the 
full SUSY algebra, e.g. by combining the Poincare and flavour symmetry group, so-called twisted 
SUSY (cf. [2] and references therein). This approach is applicable to systems with extended SUSY 
and leads to so-called 2-exact discretisations. 

Yet another difficulty, and maybe the most severe for supersymmetry on the lattice, is the 
fact that supersymmetric models with broken supersymmetry inherently suffer from a fermion sign 
problem that hinders Monte Carlo simulations of such models on the lattice. This can easily be 
seen as follows. The vanishing of the Witten index 



where F is the fermion number and H the Hamiltonian of the system, provides a necessary (but 
not sufficient) condition for spontaneous supersymmetry breaking. On the other hand, the index is 
equivalent to the partition function with periodic boundary conditions. 



and the only way for the path integral to vanish is through the fermionic determinant (or Pfaffian) 
being indefinite, independent of the fermion discretisation. Indeed, this has been seen in many 



W = lim Tr(-l)^exp(-j3//) 
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Studies of supersymmetric models on the lattice that allow spontaneous supersymmetry breaking, 
e.g. SUSY quantum mechanics with a supersymmetry breaking superpotential^ [3, 4, 5], ^ = 16 
Yang-Mills quantum mechanics [6, 7, 8, 9], ■J/' = 1 Wess-Zumino model in 2D [10], or ^ = (2,2) 
Super- Yang-Mills in 2D [11] (see, however, also [12, 13]). 

Here we propose a novel approach based on [14, 15] that circumvents the fermionic sign 
problem by formulating the path integral on the lattice in terms of fermion loops. The formulation is 
based on the exact hopping expansion of the fermionic action and allows the explicit decomposition 
of the partition function into bosonic and fermionic contributions. Consequently, one can then 
devise a simulation algorithm that separately samples the fermionic and bosonic sectors, as well as 
the relative probabilities between them. This then allows a precise calculation of the Witten index 
and a direct determination of the presence or absence of a Goldstino mode. Furthermore, although 
this is less relevant in the present context, the approach eliminates critical slowing down and also 
allows simulations directly in the massless limit or at negative bare mass values [14]. 

2. Fermion sign problem from spontaneous SUSY breaking 

Let us briefly elaborate further on the issue of spontaneous SUSY breaking (SSB), the vanish- 
ing of the Witten index and the connection to the fermion sign problem. 



It is well known that the Witten index provides a necessary, but not sufficient, condition for 
SSB [16]. One has 



From the definition it is clear that the index counts the difference between the number of hb bosonic 
and np fermionic zero energy states: 



In a field theoretic language the index is equivalent to the partition function of the system with 
periodic boundary conditions imposed on both the bosonic and fermionic degrees of freedom. 



Here the determinant has been obtained by integrating out the complex valued Dirac fermion 
fields^, while Sb{<^) is the action for the bosonic degrees of freedom, collectively denoted by ^. It 
is now clear- that in order to obtain a vanishing Witten index, we need both positive and negative 
contributions to the path integral, and this can only be achieved by the fermion determinant be- 
ing indefinite. This is the source of the fermion sign problem in the context of spontaneous SUSY 
breaking, and we argue that such a sign problem must occur in any model aspiring to accommodate 
spontaneous SUSY breaking. 

'Note that in the context of SUSY quantum mechanics it is misleading to speak of spontaneous or dynamical SUSY 
breaking; it is rather a static breaking determined by the form of the superpotential. 

^In case one is dealing with real-valued Majorana fermion fields, one obtains the Pfaffian Pf [^(0)] instead of the 
determinant. 




W= lim [TrBexp(-jS//)-Trf exp(-j8//)] 
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It is instructive to illustrate the argument in the explicit example of SUSY quantum mechanics. 
The continuum action of ^ = 2 supersymmetric quantum mechanics reads 

s=jd^\ [^)'+\p'm)f+m wit), (2.1) 

where the real field ^ denotes the bosonic coordinate, while and i/A denote the two fermionic co- 
ordinates. P{<^) is the superpotential and the derivatives P' and P" are taken with respect to 0. The 
(regulated) fermion determinant with periodic boundary conditions can be calculated analytically 
[17, 18], 

'd,+P"{<^)- 



det 



df +m 



P"(d)) 
sinh / — i^dt. 







p 

and by rewriting the sinh-function in terms of two exponentials, one can separate the positive and 
negative, or rather, the bosonic and fermionic contributions to the partition function. 



det 



dt+P"{(l>) 
df +m 



-exp + / —^dt - -exp - / —^dt =^ Zq -Zi . 



As an example, consider the superpotential Pe{(l>) = jnicj)^ + ^gcj)"^, which is even under the parity 
transformation — = — 0. In this case one finds Pg{<p) > and hence Zq / Zi, i.e. no SUSY 
breaking. On the other hand, for the superpotential /'o(0) = + ^^<P^> which is odd under 
parity, one has P''{^) = —P''{<p) and, since Sb{^) = Sb{<I>), one finds Zq = Zi, i.e. a vanishing 
Witten index and the corresponding SUSY breaking. Here, the vanishing of the Witten index is 
guaranteed by the fact that to each configuration there exists another configuration (the parity 
transformed one) that contributes to the path integral with the same weight, but with opposite sign 
stemming from the fermion determinant. Furthermore, Zq = Zi means that (in the limit of zero 
temperature, i.e. j3 — oo) the free energies of the bosonic and fermionic vacua are equal, and that 
proves the existence of a massless, fermionic mode connecting the two vacua, i.e. the Goldstino 
mode. 

Turning now to SUSY quantum mechanics on the lattice one obtains with a Wilson type dis- 
cretisation (cf. next section for further details) 

det [V* +P"(0)]p = n [1 +P"i^')] - 1 > (2-2) 
t 

where t now denotes a discrete lattice site index and V* is the backward derivative. Also in this 
case one can identify the bosonic and fermionic contributions (i.e. the first and second term of the 
difference in eq.(2.2)), and we will show below that this separation is always explicit in the fermion 
loop formulation. As a side remark, let us note that in the limit of zero lattice spacing one finds 

lini det [V* +P"] exp (^+ J^^ ^^^dt^ det [dt +P"((^>)] , 

where the exponential term can be understood as coming from radiative contributions that need 
to be corrected by 'fine-tuning' a corresponding counterterm [19, 18]. Reconsidering the two 
examples for the superpotential mentioned above, we find for Pg with m > and g > that 

det [V* + P'J] = Y{[\+m + 3g^f] - 1 > , 
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while for Pa one finds 

det [V* = n [1 - 1 , (2.3) 

t 

which turns out to be indefinite, even when A > 0. While this is necessary in order to enable a 
vanishing Witten index, it imposes a serious problem on any Monte Carlo simulation, for which 
positive weights are strictly required. Moreover, the sign problem is severe in the sense that towards 
the continuum limit (i.e. when the lattice volume goes to infinity), the fluctuations of the first 
summand in eq.(2.3) around 1 tend to zero, such that — )• is exactly realised in that limit. Hence, 
the source of the fermionic sign problem lies in the exact cancellation of the first and the second 
summand in eq.(2.3), i.e. of the bosonic and fermionic contribution to the partition function, and 
this observation also holds more generally in higher dimensions. 

In the loop formulation, to be discussed in the next section, the separation of the partition 
function into the various fermionic and bosonic sectors is made explicit and allows the construction 
of a simulation algorithm that samples these sectors sepai^ately, and more importantly also samples 
the relative weights between them. In this way, the loop formulation eventually provides a solution 
to the fermion sign problem. 

3. Loop formulation and separation of fermionic and bosonic sectors 

In this section we illustrate the loop formulation and the separation of the partition function 
into its fermionic and bosonic sectors by means of the = 1 Wess-Zumino model in two dimen- 
sions and = 2 supersymmetric quantum mechanics in one dimension. 

3.1 = \ Wess-Zumino model in two dimensions 

The Lagrangian of the .y^ = 1 Wess-Zumino model in two dimensions is given by 

^=\[d^^f + \p'{<^f + \^f[^ +P'\^))V, (3.1) 

where i// is a real, two-component Majorana field, <p a real bosonic field and P{(p) an arbiti'aiy 
superpotential. Integrating out the fermionic Majorana fields yields a Pfaffian which in general, as 
discussed above, is not positive definite. 

On the lattice, one can use the exact reformulation of the fermionic Majorana degrees of free- 
dom in terms of non-intersecting, self-avoiding loops, in order to separate the contributions of the 
Pfaffian to the various bosonic and fermionic sectors of the partition function. A similar exact re- 
formulation of the bosonic degrees of freedom in terms of bonds, can also be accomplished [20]. 
While this is not necessary for the solution of the sign problem, it provides a convenient way to 
simulate also those degrees of freedom without critical slowing down, and hence we will discuss 
this construction below. 

Employing the Wilson lattice discretisation for the fermionic part of the Lagrangian in eq.(3. 1) 
yields 
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where t, now represents the real, 2-component Grassmann field, while 'lo = —'lo^ is the charge 
conjugation matrix and V* the Wilson term. Using the nilpotency of Grassmann elements one 
can expand the Boltzmann factor leading to 



n ( 1 - \m{<^.)^i^£,. ) n(i + 



where M{<^x) = 2 + P"(0,v)> r(ib/x) = ^(1 =F 7^j) and x denotes the discrete lattice site index. Per- 
forming now the integration over the fermion field, at each site x the fields <^J^ and t,x must be 
exactly paired in order to give a contribution to the path integral, so one finds 

where the occupation numbers m{x) =0,1 for the monomer terms and (x) =0,1 for the fermionic 
bonds (or dimers), satisfy the constraint 

m{x) + \Y,{b^{x)+b^{x-fi)) = \ Wx. (3.2) 

This constraint is equivalent to the fact that only closed, self-avoiding paths survive the Grass- 
mann integration. When integrating out the fermion fields, the projectors r(jLt) eventually yield 
a weight CO, which only depends on the geometric structure of the specific constrained path (CP) 
configuration £ G In particular, one has 

|a)(^)| =2-"^/2, 

where Uc denotes the number of corners in the loop configuration, while the sign depends on the 
topology of the loop configuration and will be discussed below. 

As mentioned above, the bosonic fields can be treated analogously [20]. On the lattice, the 
kinetic term {d^<p)^ yields ^x'Px-fi^ and expanding this hopping term to all orders gives 

^•^n I (•^.v^.-a)"'^'^ n^'^p (-^^vi^^)) M^'"^'^ (3-3) 

with bosonic bond occupation numbers n^{x) = 0, 1 , 2, . . . and V (0^) = 4 + P'{<^x)^. In contrast to 
the fermionic case, the exact reformulation requires one to include an infinite number of terms in 
the hopping expansion, and hence occupation numbers up to infinity, instead of just and 1 as for 
the fermionic bonds. Integrating out the bosonic fields yields the site weights 



e(n^(x),m(x)) = j d^, exp (-^-V{^,)\ ^f'^M{^, 



where N{x) = {x) +n^{x — jX)) counts the number of bosonic bonds attached to a given site, 

while M((^;f)'"(^) may contribute additional powers of So the bosonic contribution to the weight 
of a given configuration factorises into a product of local weights, 

W{n^{x)Mx)) = n ^TTne('V • 
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In summary, the fermionic and bosonic degrees of freedom in our original partition function 
are now expressed in terms of fermionic monomers and dimers and bosonic bonds, and the integra- 
tion over the fields has been replaced by a constrained sum over all allowed monomer-dimer-bond 
configurations, yielding 

Z= E E \W{n^{x)M^))-m\- 

In particular, the partition function for the fermionic degrees of freedom is represented by a 
sum over all non-oriented, self-avoiding fermion loops 

E \W{n^{x),m{x))-(0i£)l if G ifooUifioUifoi Uifu , (3.4) 

where i represents a fermion loop configuration in one of the four topological classes , with 
h,h = 0,1 denoting the total number of loop windings (modulo 2) along the first and second 
direction, respectively, on the periodic lattice torus. CO denotes the weight of the specific loop con- 
figuration and depends on the geometry of the loop configuration, as discussed above. The sign 
of the weight is solely determined by the topological class and the fermionic boundary conditions 
e^, where £^ = 0, 1 stands for periodic and anti-periodic boundary conditions, respectively. Con- 
figurations in have positive weights independent of the boundary conditions, while the sign of 
the weights for configurations in Jfoii^io and is given by (— l)''' '^f So the partition func- 
tion in eq.(3.4), where all sectors contribute positively, represents a system with unspecified 
fermionic boundary conditions [21], while a partition function with fermionic b.c. periodic in the 
spatial and anti-periodic in the temporal direction, respectively, is described by the combination 

This combination represents the system at finite temperature. Analogously, the partition function 
with fermionic b.c. periodic in all directions - the Witten index - is given by 

The interpretation of the Witten index in terms of the partition functions Zc^.. is straightforward. 
Any fermion loop winding non-trivially around the lattice carries fermion number F = \, hence 
configurations with an odd number of windings, i.e. configurations in Z_^^^,Zj^^^ and Z_sfu, also 
carry fermion number F = \, while configurations with no, or an even number of windings, i.e. in 
Z_5fQQ, have F = 0. The partition function Z^^^^ may therefore be interpreted as representing the 
bosonic vacuum, while the combination Z^,q +Z_^q| +Z^^^ corresponds to the fermionic vacuum. 
Consequently, the latter contributes to the Witten index W = Zpp with opposite sign relative to the 
bosonic vacuum. Since each of the four partition functions is positive, vanishing of the Witten 
index implies Z^^ = Z^,o + Z^,^ 

3.2 = 2 supersymmetric quantum mechanics 

The loop formulation for the .yV = 2 supersymmetric quantum mechanics on the lattice is 
obtained in a completely analogous manner. Using again the Wilson lattice discretisation for the 
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fermionic part, the continuum action in eq.(2.1) reads 




where x now denotes the one-dimensional, discrete lattice site index. Note that in one dimension the 
fermionic lattice derivative, including the contribution from the Wilson term with Wilson parameter 
r = 1, becomes a simple, directed hop YxWx-i which, in the loop formulation, can be described 
by the (directed) bond occupation number b{x) = 0, 1. Integrating out the fermionic degrees of 
freedom yields a constraint for the fermion monomer and bond occupation numbers, analogous to 
eq.(3.2), namely 

m{x) + ^{b{x)+b{x-l)) = l \/x. 

So for .A' = 2 supersymmetric quantum mechanics the loop formulation becomes particularly sim- 
ple: there are just two different loop configurations, namely a bosonic one, where m{x) = 1 ,b{x) = 
for all X, and a fermionic one where m{x) = 0,b{x) = 1 for all x. Since the latter corresponds to a 
closed fermion loop, it will pick up a minus sign from the Grassmann integration, relative to the 
bosonic contribution. 

On top of the two fermion loop configurations one may treat the bosonic fields in the same 
way as before and employ a hopping expansion to all orders. After rearranging the bosonic fields 
the integration can eventually be performed separately at each site and one ends up with the weight 

W{n{x)M^))=Y{^, I d4x ^^^^''^^-'^ e-J^*-)(l+P"(,/.,))'"« (3.5) 
\^n[x)\ J 

for a given bosonic and fermionic bond configuration, with V{(^x) = 2 + P'{<p)^. In terms of these 
weights the partition function can now be written as 

Z^= f ^(p^xj/^Xj/Q-^^ = L \W{n{x),m{x))\, ^e^oU^i, (3.6) 

where the second sum is over all allowed bosonic bond configurations {n} G CP, and £ represents 
one of the two fermion loop configurations in the topological classes or , respectively. As 
before, the sign of the weight is solely determined by the topological class and the fermionic 
boundary condition. If / = 0, 1 denotes the fermion loop winding number and, as before, e = 0, 1 
the periodic and anti-periodic fermionic boundary condition, respectively, the sign is given by 
(_l)'-(e+l). 

Choosing anti-periodic fermionic boundary conditions £ = 1 we find 

which is simply the partition function for the system at finite temperature, while choosing periodic 
fermionic boundary conditions £ = yields 

representing the Witten index. Here, the interpretation is particularly intuitive: the two fermion 
loop configurations simply represent the bosonic and fermionic vacuum, while Z and Z repre- 
sent the bosonic and fermionic partition function in the corresponding sectors. The Witten index 
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vanishes whenever =Z^^, i.e. when the contributions from the bosonic and fermionic sectors 
cancel. In this case, the free energy of the bosonic and fermionic vacuum is equal, and this is 
equivalent to saying that there exists a gapless, fermionic excitation which oscillates between the 
two vacua, i.e. the Goldstino mode. 

4. Solution of the fermion sign problem 

In this section we briefly describe the simulation algorithm and explain how it eventually 
solves the fermion sign problem. The loop system can most efficiently be simulated by enlarging 
the configuration space by open strings. Following [20] the bosonic bonds are updated by inserting 
two bosonic sources, which sample directly the bosonic 2-point correlation function. Similarly, 
the fermion bonds are most efficiently updated by simulating a fermionic string [14] that sam- 
ples the configuration space of the fermionic 2-point correlation function, instead of the standard 
configuration space of loops. 

For the o/K = 1 Wess-Zumino model in two dimensions, where the fermionic degrees of free- 
dom are Majorana, the open fermionic string is non-oriented and corresponds to the insertion of a 
Majorana fermion pair {i^J*^, ^y} at position x and y, while for ^ = 2 supersymmetric quantum 
mechanics, where the fermionic degrees of freedom are Dirac, the string is oriented and corre- 
sponds to the insertion of a Dirac fermion pair {i//^, 

The algorithm proceeds by locally updating the endpoints of the open fermionic string using a 
simple Metropolis or heat bath step according to the weights of the corresponding 2-point function 
(cf. [14] for details). When one end is shifted from, say, x to one of its neighbouring sites y, a 
fermionic dimer on the corresponding bond is destroyed or created depending on whether the bond 
is occupied or not. Contact with the partition functions Zy.. and Z^., respectively, is made each 
time the open string closes. This then provides the proper normalisation for the expectation value 
of the 2-point function. 

The solution of the fermion sign problem discussed above relies on the correct determination 
of the relative weights between the bosonic and fermionic sectors, and the simulation algorithm 
described in [14] achieves this in a most efficient way. The open fermionic string tunnels between 
loop configurations in the various topological homotopy classes ^oo^-^^'io, ^01,^11 in two, and 
in one dimension, thereby determining the relative weights between the partition functions 
Z , Z , Z_sfQj , Z_^, , , or Z_^, , Z_^, , respectively. From the relative weights, the Witten index (or 
any other partition function of interest) can be reconstructed a posteriori. 

Let us emphasise that the open string algorithm and the corresponding absence of critical slow- 
ing down at the critical point as reported in [14] is crucial for the solution of the sign problem. Since 
the algorithm updates the configurations according to the fermionic 2-point correlation function, 
they are updated equally well on all length scales up to a scale of the order of the largest fermionic 
correlation length. This is in fact the reason why critical slowing down is essentially absent even 
at a critical point when the correlation length becomes infinite. Now, in order to have IV = in the 
continuum (L — )• 00) the Goldstino mode has to become massless. Since the algorithm ensures an 
efficient update of that mode, the tunneling between the bosonic and fermionic vacua is guaranteed 
and the Witten index indeed vanishes in practice. 
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Figure 1: Spectrum for the perturbatively improved standard discretisation at g/m'^ ~ 1.0. In the left plot 
the excitation energies mg (circles) and mf (squares) and the lattice spacing are expressed in units of the 
lattice extent L, while in the right plot they are expressed in units of the bare mass parameter m. 



5. Results ioY ,yV = 2 supersymmetric quantum mechanics 

We are now ready to present some results for the spectrum and the Witten index for the case 
of .yV = 2 supersymmetric quantum mechanics. 

5.1 Spectrum for the perturbatively improved standard discretisation 

Here we present our results for the spectrum using the standard discretisation including a 
counterterm. As discussed in [19], for the standard discretisation described above, the coiTcct 
continuum limit is spoiled by radiative corrections^. This can be corrected by adding a counter- 
term of the form jLv^'l^-v) [19]. By doing so, one ensures that all observables reach the coiTcct 
continuum limit and that the full supersymmetry is eventually restored. 

As an example we consider a superpotential with unbroken supersymmetry, i.e. Pe{^) = 
^m(l)^ + and choose the coupling g/rri^ = 1.0. The results are presented in figure 1 where 

we show the lowest lying excitation energies for the boson (circles) and the fermion (squares) as 
a function of the lattice spacing a for various values of fixed mL. In the left plot, the quantities 
are expressed in units of the lattice extent L, while in the right plot, they are expressed in units of 
the bare mass parameter m in order to illustrate the common scaling behaviour. The leading lattice 
artifacts turn out to be i^{a) for both the fermion and boson masses. At finite lattice spacing the 
supersymmetry is explicitly broken by the discretisation, and hence the boson and fermion masses 
are not degenerate. In the continuum limit, however, the supersymmetry is restored and the masses 
become degenerate. 

5.2 Spectrum for the 2-exact discretisation 

As briefly discussed in the introduction, for models with extended supersymmetry it is some- 
times possible to preserve some of the supersymmetries exactly at finite lattice spacing [22] . The 
so-called Q-exact discretisations preserve a suitable sub-algebra of the full supersymmetry algebra, 

^Note, however, that in one dimension these corrections are finite. 
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Figure 2: Spectram for the Q-exact discretisation at g/m'^ = 1.0. In the left plot the excitation energies mg 
(circles) and mf (squares) and the lattice spacing are expressed in units of the lattice extent L, while in the 
right plot they are expressed in units of the bare mass parameter m. The dashed line denotes the exact value 
obtained with Numerov's method. 



1. e. a linear combination of the available supersymmetries. In the context of ,yV = 2 supersymmet- 
ric quantum mechanics the 2-exact action is obtained from the standard action by adding, e.g., the 
term Y.xP' {^x)^* <^x, but other forms are also possible [23, 18]. Since the term contains a derivative, 
there are additional hopping terms that need to be considered in the hopping expansion [24]. 

In the following we concentrate again on the superpotential with unbroken supersymmetry. 
Using the 2-exact discretisation one expects degenerate fermion and boson masses even at finite 
lattice spacing [22] and this is beautifully confirmed by our results at g/m^ = 1 .0 presented in figure 

2. Note that the leading lattice artifacts are again ff(a) for both the fermion and boson masses. 

Finally, in the right plot we show the results of a high precision simulation that serves the 
purpose of checking the coiTcctness of the new simulation algorithm, as well as our procedures for 
the extraction of the fermion and boson masses. Indeed, we can confirm the mass degeneracy to 
a precision better than a few per mill at all lattice spacings, and the continuum value of the mass 
gap agrees with the exact result in the continuum obtained with Numerov's method (dashed line) 
also within a few per mill. Due to the loop formulation and the efficiency of the new simulation 
algorithm, these results can be obtained with a very modest computational effort. 

5.3 Witten index 

Let us turn to the Witten index W °< Zp/Z^, i.e. the partition function with periodic boundary 
conditions Zp, normalised to the finite temperature partition function Z^. We start with the super- 
potential Pe{<^) = \m<^^ + \g<^'^ for which supersymmetry is unbroken and W ^0. The results ai^e 
presented in figure 3, where we show Zp/Za in the left plot as a function of the bai^e mass parameter 
am for the coupling g/rn^ = \.0 using the Q-exact discretisation. The continuum limit is reached as 
am — >• 0, and we indeed find that Zp/Za — ^ 1> i-C- IV / in that limit on a sufficiently large lattice. 
The fact that at fixed lattice extent L the partition function ratio goes to zero with am — )• can be 
interpreted as a 'finite size', or rather finite temperature effect, since the temperature T is inversely 
proportional to the extent of the lattice L. So if we plot the data as a function of m/T = mL, mL = 
corresponds to infinitely high temperature, while the limit mL — )• oo corresponds to T — )• 0. This is 
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am mL 



Figure 3: Simulation results for the Witten index W °^ Zp/Za as a function of the bare mass parameter 
am (left plot), and as a function of the inverse temperature mL (right plot), for a system with unbroken 
supersymmetry using the g-exact discretisation at g/m^ — 1.0. 




Figure 4: Simulation results for the Witten index W °^ Zp/Z^ as a function of the inverse temperature mL 
(left plot), and as a function of the lattice spacing a /L at fixed temperature (right plot), for a system with 
broken supersymmetry using the standard, perturbatively improved discretisation at A^/ot^ — 1.0. 



illustrated in the right plot of figure 3, from where we find Zp/Za — )■ 1 in the continuum limit for 
temperatures mL > 5. The data shows a rather good scaling behaviour towards the continuum limit 
(L — )• oo at fixed mL). In the continuum, the value Zp/Za = 1 is approached exponentially fast with 
mL, i.e. towards zero temperature mL oo. 

Note that with the formulation and algorithm presented here, it is also possible to simulate at 
negative bare mass m < 0. In this case, we still find Zp/Za — 1 in the continuum limit towards zero 
temperature, albeit at a slower rate. 

When repeating this exercise for the superpotential Po(</>) = —^<P + ^^'P^> for which su- 
persymmetry is broken, we should expect a vanishing Witten index. Our results for this case are 
presented in figure 4, where we show the Witten index W oc Zp/Za as a function of mL in the left 
plot, and as a function of the lattice spacing a/L at fixed values of mL in the right plot using the 
standard, perturbatively improved discretisation at X^/m^ = 1.0. In this case we find that while the 
Witten index W Zp/Za approaches zero at 'infinite temperature' mL — )• as in the unbroken case, 
it also does so at any value of mL in the continuum limit, even though for large values of mL the 
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scaling towards the continuum limit is reached only at very fine lattice spacings. In fact, for large 
mL the approach W — )• is exponentially slow in the lattice spacing. This is exemplified in the 
right plot of figure 4, where the dashed lines are the results of an analytic calculation which will be 
reported elsewhere [24]. The absence of critical slowing down for our update algorithm guarantees 
that reliable results can be obtained despite the exponentially slow approach to the continuum. 

6. Summary and outlook 

We have discussed the occurrence of a fermion sign problem in the context of spontaneous 
supersymmetry breaking on the lattice and its relevance for the vanishing of the Witten index, 
regulated as a path integral on the lattice. We then argued that with the help of the fermion loop 
expansion one can achieve an explicit separation of the bosonic and fermionic contributions to the 
path integral in such a way that the source of the sign problem, namely the cancellation between the 
bosonic and fermionic contributions to the partition function with periodic boundary conditions, is 
isolated. The solution of the fermion sign problem is then achieved by devising an algorithm 
which separately samples the bosonic and fermionic contributions to the partition functions and, in 
addition, also samples the relative weights between them, essentially without any critical slowing 
down. In such a way one is able to calculate the Witten index on the lattice without suffering from 
the fermion sign problem even when the index vanishes. The absence of critical slowing down is 
essentially due to the fact that the algorithm directly samples the massless Goldstino mode which 
mediates the tunnelings between the bosonic and fermionic vacua. 

As examples we described in some detail the exact reformulation of the lattice path integral 
in terms of fermionic bonds and monomers, and bosonic bonds for the o/K = 1 supersymmetric 
Wess-Zumino model in 2 dimensions and .yV = 2 supersymmetric quantum mechanics. For the 
latter we presented some results from lattice simulations. For a superpotential with unbroken su- 
persymmetry, we calculate the energy gap of the lowest bosonic and fermionic excitation using 
both the standard discretisation with a fine-tuned counterterm, as well as a Q-exact discretisation 
which preserves a linear combination of the two supersymme tries. In the latter case the boson and 
fermion spectra ai^e degenerate even at finite lattice spacing, while in the former case they become 
degenerate only in the continuum limit. Finally, we also present lattice calculations of the Witten 
index for broken and unbroken supersymmetry using the standard discretisation with a counter- 
term. For both cases we are able to reproduce the correct Witten index in the continuum limit. 
For broken supersymmetry the approach to the continuum limit is exponentially slow in the lattice 
spacing. Using the loop formulation and the fermion worm algorithm [14] the exponentially slow 
approach as well as the sign problem is no obstacle in practice. 

Obviously, the approach presented here is particularly interesting for the -yV = 1 supersym- 
metric Wess-Zumino model in 2 dimensions, where so far simulations on the lattice have suffered 
from the fermion sign problem [10]. Work in this direction is in progress. 
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